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Laminar and turbulent perfect-gas air, Navier-Stokes computations have been performed for a pro- 
posed Mars Smart Lander entry vehicle at Mach 6 over a free stream Reynolds number range of 6.9xl0 6 /m to 
2.4x1 0 7 /m (2.1xl0 6 /ft to 7.3xl0 6 /ft) for angles-of-attack of 0-deg, 11-deg, 16-deg, and 20-deg, and comparisons 
were made to wind tunnel heating data obtained at the same conditions. Boundary layer edge properties were 
extracted from the solutions and used to correlate experimental data on the effects of heat-shield penetrations 
(bolt-holes where the entry vehicle would be attached to the propulsion module during transit to Mars) on 
boundary-layer transition. A non-equilibrium Martian-atmosphere computation was performed for the peak 
heating point on the entry trajectory in order to determine if the penetrations would produce boundary-layer 
transition by using this correlation. 


Nomenclature 

a w speed of sound at wall (m/s) 

D diameter (m) 

H aw adiabatic enthalpy (J/kg) 

H w wall enthalpy (J/kg) 

H 0 total enthalpy (J/kg) 

li heat transfer coefficient (kg/m 2 /sec). 

h = <■//(//.,,. ■//.,) 

hpp Fay-Riddell heat transfer 

coefficient (kg/m 2 /sec) 

M„ free stream Mach number 

M e boundary layer edge Mach number 

free stream pressure (N/tn 2 ) 
q heat transfer rate (W/rn 2 ) 

r radial position (m) 

R base (maximum) radius (m) 

R„ nose radius (m) 

R c comer radius (m) 

free stream unit Reynolds number (I/m) 
(Re we ) mrh critical Reynolds number for turbulence, 
(*e«.e) tllrh = 

Re ce n wall cell Reynolds number, 

Re ceU = (Pk'V^'X 

Reg boundary layer momentum thickness 

Reynolds number, Re e = (p l ,U t ,6)/\i jl . e 
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T „ free stream temperature (K) 

U„ free stream velocity (m/sec) 

M’ penetration diameter (cm) 

a angle-of-attack (deg) 

8 boundary layer thickness (m) 

ip model rotation (deg) 

As wall cell height (m) 

wall cell viscosity (kg/m/s) 
p„ free stream density (kg/m 1 ) 

wall cell density (kg/m 1 ) 

Introduction 

The Mars Smart Lander (MSL) mission (Refs. 
1 and 2) is intended to demonstrate an entry into the 
Martian atmosphere with an uncertainty in landing coor- 
dinates of less than 10 km. One of the important design 
issues for this vehicle is the determination of whether 
the boundary layer on the forebody of the entry vehicle 
will be laminar or turbulent. This determination is more 
complicated than has been the case for previous Mars 
missions because the MSL forebody will not have a 
smooth, unbroken surface; instead there will be six, cir- 
cular penetrations in the forebody heat shield where 
bolts will be used to attach the vehicle to a cruise stage 
during transit from Earth. 

The goal of the present study was to formulate 
transition criteria for the MSL which account for pene- 
tration diameter, penetration location, angle-of-attack, 
and Reynolds number. To accomplish this goal, wind 
tunnel testing was conducted in order to generate a tran- 
sition database for these parametrics, and then 
Navier-Stokes computations were performed for the 
wind tunnel conditions in order to determine bound- 
ary-layer edge quantities which were used to correlate 
the transition data. The computations and correlations 


1 

American Institute of Aeronautics and Astronautics 



are presented in this paper, while the wind tunnel data 
are presented in the companion paper by Liechty in Ref. 
3. Comparisons of both laminar and turbulent aeroheat- 
ing predictions with the experimental data are also pre- 
sented, as are predictions for the peak heating condition 
during Mars entry. 

A detailed investigation of the MSL aeroheat- 
ing environment during Mars entry is presented by 
Edquist in Ref. 4. A similar study on penetration effects 
was performed for the Cienesis mission, and is reported 
in Ref. 5. 





Figure 1: Dimensions of Full-Scale Mars Smart 
Lander Entry Vehicle 



Figure 2: Mars Smart Lander and Cruise Stage 


Mars Smart Lander Geometry 

The proposed Mars Smart Lander entry vehicle 
is a 70-deg half-angle sphere-cone with a biconic after- 
body (Fig. 1). The forebody will have six penetrations 
spaced at 60-deg increments where it will be bolted to 
the carrier vehicle in transit to Mars (Fig. 2). In this 
study, various penetration diameters and radial and 
angular positions were investigated in order to deter- 
mine their effects on boundary layer transition and sur- 
face aeroheadng levels. 

Computational Method 

Computations were performed using the 
LAURA (Refs. 6 and 7) code (version 4.9.2). The 
LAURA (Langley Aerothermodynamic Upwind Relax- 
ation Algorithm) code is a three-dimensional, finite-vol- 
ume solver which includes perfect-gas, equilibrium and 
non-equilibrium chemistry models. The code can be 
used to solve the inviscid, thin-layer Navier-Stokes, or 
full Navier-Stokes equations. For the current study the 
thin-layer mode was employed; it was concluded in Ref. 
5 from computations on a similar blunt body that this 
mode provided accurate results for attached forebody 
Hows. Time integration to steady-state in LAURA is 
accomplished through a point-relaxation scheme. 
Roe-averaging (Ref. 8) with Flarten’s entropy fix (Ref. 
9) and Yee’s Symmetric Total Variation Diminishing 
limiter (Ref. 10) is used for inviscid fluxes, and a sec- 
ond-order scheme is employed for viscous fluxes. For 
turbulent computations, the algebraic Baldwin-Lomax 
(Ref. 11) model with modifications (Ref. 12) for com- 
pressible flow and the Dhawan-Narashima (Ref. 13) 
transition model were employed. In this study, the per- 
fect-gas air model was used for the wind tunnel compu- 
tations, and an 8-species non-equilibrium, non-ionizing 
Martian atmosphere model (C0 2 , CO, N 2 , 0 2 , NO, C. 
N, O) was used for the flight case. 

Free stream conditions for the LAURA wind 
tunnel computations were set to the nominal free stream 
operating conditions of the NASA Langley Research 
Center (LaRC) 20-Inch Mach Air Tunnel, which are 
listed in Table 1. For the wind tunnel computations, a 
uniform, ambient 300 K wall temperature boundary 
condition was imposed. The use of a constant wall tem- 
perature is valid because the experimental data are 
reported in terms of the non-dimensional ratio, h/h FR , 
which is assumed to remain constant with wall tempera- 
ture. The quantity h FR is the heat-transfer coefficient 
computed using the Fay-Riddell (Ref. 14) method with 
the same nose radius as the MSL model at a wall tem- 
perature of 300 K (340 °R). 
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For the flight case, free stream conditions were 
taken from the peak heating point on the nominal trajec- 
tory and are listed in Table 2. The wall boundary condi- 
tion was set to “super-catalytic” (full recombination to 
free stream concentrations) with a radiative equilibrium 
wall temperature and a wall emissivity of 0.90. 

A structured, finite-volume, multiple-block 
grid with a singularity-free nose was employed for the 
computations. Although the grid completely encom- 
passed the afterbody and wake, computations were lim- 
ited to the forebody blocks as this was the region of 
interest in the present study. The forebody grid blocks 
contained approximately 230,000 points with a 
body-normal (/.--index) resolution of 65 points. Grid 
adaptation was performed (as per the method detailed in 
Ref. 7) to align the grid with the bow shock and to pro- 
duce nominal wall cell Reynolds numbers on the order 
of Re cell = 10. 

The effects of normal grid-point resolution on 
the computed heating distributions was examined by 
repeating the computations for the a = 16-deg, Re x = 

1.9xl0 7 /m case with grids containing half (32) and twice 
(128) as many cells in the normal direction as the nomi- 
nal grid (64). The heating distributions from these com- 
putations are shown in Fig. 3. Heating levels dropped 
nearly uniformly over the entire surface by about 5% 
from the 32-cell grid to the 64-cell grid, but from the 
64-cell grid to the 128-cell grid, heating levels dropped 
by about 1% except around the stagnadon region, where 
the decrease was approximately 2%. Therefore, it was 
concluded that the original 64 normal cell grid provided 
acceptable accuracy for this study. 

Although the wind tunnel test models were fab- 
ricated with actual penetrations to replicate those in the 


flight vehicle heat shield, the MSL computational grid 
did not include these penetrations. Thus, the effects of 
the penetrations on the flow field were not modeled in 
the computations. The effects of the penetrations on the 
state of the boundary layer were simulated by specifying 
that transition began at the penetradon location. In 
order to attempt to bound the effects of the penetrations 
on boundary layer transition, the transition length in the 
Dhawan-Narashima transition model was specified as 
either zero, to simulate an immediate jump to 
fully-developed turbulent flow, or as equal to the run- 
ning length of the boundary layer from the nose of the 
vehicle to the penetration location to simulate natural 
transition. 



1.00 -0.50 0.00 0.50 1.00 


r/H 

Figure 3: Effects of Normal Grid Resolution on 
Computed Heating Levels, a = 16, Re„ = 2.4xl0 7 /m 


Table 1: Free Stream Conditions for Wind Ttinnel Cases 


Re x 

(1/m) 

Moo 

Too 

(K) 

Poo 

(kg/m 3 ) 

Uoo 

(m/s) 

h FR 

(kg/m 2 -s) 

Qfr 

(W/cm 2 ) 

6.7xl0 6 

5.95 

62.0 

3.35x1 O' 2 

938.6 

0.283 

5.70 

8.5xl0 6 

5.97 

62.2 

4.05xl0" 2 

943.0 

0.313 

6.45 

9.8xl0 6 

5.98 

62.2 

4.62xl0" 2 

944.4 

0.335 

6.95 

l.lxlO 7 

5.99 

61.6 

5.29xl0" 2 

940.1 

0.356 

7.25 

1.4xl0 7 

6.00 

61.3 

6.41xl0" 2 

940.4 

0.392 

7.95 

1.7xl() 7 

6.02 

63.4 

7.92xl0" 2 

958.7 

0.446 

9.92 

1.9xl0 7 

6.03 

62.8 

8.99x1 0 -2 

955.6 

0.474 

10.4 

2.4xl0 7 

6.06 

62.3 

1. 13x10-* 

954.6 

0.529 

11.5 
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Table 2: Free Stream Conditions at Peak Heating 


Time (sec) 

163 

Alt. (km) 

37.1 

(m/s) 

4919 

p«,(kg/m 3 ) 

4.15x1 0 -4 

Too (A) 

159.3 

Moo 

24.7 

Rejl/m) 

2.57 x 10 4 5 



Outer 
Radial 
Location 
r/R = 0.70 


r/R = 


Inner 

Radial 

Location 

(Penetrations not 
shown to scale) 


Penetration Sizes 

* 7.62 cm (3.0-in.) lull scale 

full scale 

• 3.81 cm (1.5 -in.) full scale 


Figure 4: Layout of Penetrations on MSL Model 


Experimental Method 

The wind tunnel test which complements this 
computational study is presented in detail in Ref. 3. and 
a brief discussion is included herein for completeness. 

Aeroheating tests were conducted in the NASA 
LaRC 20-Inch Mach 6 Air Tunnel. This facility is a 
blow-down tunnel in which heated, dried, and filtered 
air is used as the test gas. The tunnel has a two dimen- 
sional, contoured nozzle which opens into a 0.521 m x 
0.508 m (20.5-in. x 20.0-in.) test section. The tunnel is 
equipped with a bottom-mounted injection system 
which can transfer a model from the sheltered model 
box to the tunnel centerline in less than 0.5 sec. Run 
times of up to 15 minutes are possible in this facility, 
although for the current aeroheating study, run times of 
only a few seconds were required. The nominal reser- 
voir conditions of this facility are stagnation pressures 
of 206.8 to 3447.4 kPa (30 to 500 psia) with stagnation 
temperatures of 422.2 to 555.5 K (760 "R to 1000 ”R), 


which produce perfect-gas free stream flows with Mach 
numbers between 5.8 and 6.1 and Reynolds numbers of 
1.64xl0 6 m’ 1 to 23.3x10 s in 1 (0.5xl0 6 ff 1 to 7.3xl0 6 
ft" 1 ). A more detailed description of this facility is pre- 
sented in Refs. 15 and 16. Representative flow condi- 
tions for each of the standard 20-Inch Mach 6 Air 
Tunnel operating points have been computed using the 
GASPROPS (Ref. 17) code and are listed in Table 1. 

Global surface heating distributions were 
obtained using the digital optical measurement method 
of two color, relative-intensity, phosphor thermography 
(Refs. 18-21). In this method, ceramic wind tunnel 
models are coated with a phosphor compound which 
fluoresces in two separate regions (green and red) of the 
visible light spectrum. Before and during a wind tunnel 
run. the phosphor-coated model is illuminated by ultra- 
violet (UV) light sources, and the resulting fluorescent 
intensity of the model is recorded and digitized through 
a three-color CCD (charge coupled device) camera. 
Intensity data are converted to surface temperature val- 
ues using system calibrations. Global heat-transfer dis- 
tributions are then computed from these temperature 
data using one-dimensional, constant heat-transfer coef- 
ficient, conduction theory. As discussed in Ref. 3, the 
estimated experimental uncertainty of the heating data is 
approximately ±13%. 

The heating data are presented in the 
non-dimensional form, h/hpp. where hpp is the refer- 
ence heat-transfer coefficient from Fay-Riddell theory 
for a 300 K (540 ”R) surface temperature and the nose 
radius of the model. In the definitions of both h and 
h FR , the adiabatic wall enthalpy, H aw is assumed to be 
equal to the total enthalpy, H 0 . 

Heating distributions were measured on 
0.0314-scale ceramic models of the MSL entry vehicle. 
Model parametrics were: penetration diameter (w = 7.62 
cm / 3.0-in. 5.59 cm / 2.2-in., or 3.81 cm / 1.5-in. full 
scale); angular location of the penetration (with respect 
to the leeside centerline); and radial location of the pen- 
etration ( r/R = 0.41 or 0.70). These parametrics are 
shown in Fig. 4. Note that the actual flight vehicle 
would have six penetrations of uniform diameter at a 
constant radial location separated by 60-deg rotational 
increments. 


Aeroheating Predictions and Com- 
parisons with Experimental Data 

Laminar Comparisons 

Laminar computations were performed at the 
conditions listed in Table 1 for angles-of-attack of 
0-deg, 11-deg, I6-deg, and 20-deg. Symmetry plane 
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comparisons of these computed heating distributions 
with the wind tunnel data for each angle-of-attack are 
presented in Figs. 5-8. In these figures, the experimental 
data shown were measured on models without penetra- 
tions on the symmetry line. The computed laminar heat- 
ing distributions were found to agree with the 
experimental data to within the estimated uncertainty for 
all cases except at the forebody comers. 

Turbulent Comparisons 

Turbulent computations were performed for 
wind tunnel conditions of Re ^ = 1.4d0 7 /m to 

2.4xl0 7 /m at a = 16-deg, which is the nominal 
angle-of-attack at peak heating on the flight trajectory. 
Because boundary layer transition in the experiment was 
produced by penetrations in the test models (as it would 
also occur in flight) as opposed to naturally, transition in 
the computations was modeled in several different ways 
in order to attempt to model the data. Computations 
were performed with: the boundary layer fully turbulent 
over the length of the vehicle; natural transition begin- 
ning at the location of the penetration; and zero-length 
transition to fully-turbulent flow at the penetration loca- 
tion. For the natural transition computations, the transi- 
tion length was set equal to the running length of the 
flow from the nose to the penetration. 

Symmetry-plane comparisons of these compu- 
tations with the experimental data for both radial pene- 
tration locations are shown for r/R = 0.41, Re x = 

1.4xl0 7 /m, 1.9xl0 7 /m, and 2.4xl0 7 /m in Figs. 9-11 and 
for the same Reynolds numbers with r/R = 0.70 in Figs. 
12-14. In both sets of figures, comparisons are shown 
only for the lee-side of the vehicle because the 
wind-side penetration was observed to have little or no 
effect on the heating except at the penetration itself. 



-1.00 -0 50 0.00 0.50 1 00 


r/R 

Figure 5: Laminar Centerline Heating Comparison 
for a = 0-deg 



Figure 6: Laminar Centerline Heating Comparison 
for a = 11-deg 



Figure 7: Laminar Centerline Heating Comparison 
for a = 16-deg 



Figure 8: Laminar Centerline Heating Comparison 
for a = 20-deg 
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Figure 9: Turbulent Leeside Centerline Comparison 
for a = 16, Re„ = 1.4x1 0 7 /m, r/R = 0.41 
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Figure 10: Turbulent Leeside Centerline Compari- 
son for a = 16, Re„ = 1.9xl0 7 /m, r/R = 0.41 



Figure 11: Turbulent Leeside Centerline Comparl 
son for a = 16, Re x = 2.4x1 0 7 /m, r/R = 0.41 
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Figure 12: Turbulent Leeside Centerline Compari- 
son for a = 16, Re„ = 1.4xl0 7 /m, r/R = 0.70 



Figure 13: Turbulent Leeside Centerline Compari- 
son for a = 16, Re, x = 1.9xl0 7 /m, r/R = 0.70 
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As shown in Figs. 9-14, the location at which 
transition was specified to begin, and the length of tran- 
sition which was specified had significant effects on the 
computed heating levels. For both radial locations of 
transition onset, the heating levels computed by assum- 
ing fully-turbulent flow beginning at the specified radial 
position (i.e. zero-length transition) were higher than 
those computed when the flow was treated as fully tur- 
bulent from the nose of the vehicle. This difference was 
due to the fact the boundary layer thickness at a given 
location was greater for the fully turbulent flow from the 
nose because of the longer running length over which 
the turbulent boundary layer had to grow. Similarly, the 
boundary layer growth was greater when zero-length 
transition was specified at the penetration location than 
when natural transition was specified. Hence, for the 
penetration location of r/R = 0.40. the heating levels at 
the end of natural transition were higher than those for 
zero-length transition at the same location because of 
the thinner boundary layer. However, for the r/R = 0.70 
penetration case, the boundary layer was still transi- 
tional at the shoulder of the vehicle when natural transi- 
tion was specified, so comparison with the 
zero-transition length computation was not possible. 

The comparisons between experimental and 
computations results were also very dependent on tran- 
sition location and transition length, as well as on Rey- 
nolds number and penetration size. 

Consider first the comparisons for the inner 
penetration location of r/R = 0.41 (Figs. 9-11). At the 
lowest Reynolds number of Re^ = 1.4.rl0 7 /m (Fig. 9), 
the smallest penetration had no effect on the heating and 
the intermediate penetration had only a slight effect, and 
these two data sets compared well with the laminar 
computation. Although a more significant heating 
increase was produced by the largest penetration, the 
heating levels only rose gradually over the length of the 
vehicle and never reached a “plateau" typical of fully 
developed turbulent flow. For this case, the measured 
heating levels fell between the laminar and natural tran- 
sition computations. 

For r/R = 0.41 at the intermediate Reynolds 
number of Re„ = 1.9.vl0 7 /m (Fig. 10). the data for the 
smallest penetration appeared to remain laminar and 
matched the laminar computation. The heating distri- 
butions for the two larger penetrations had shapes simi- 
lar to the natural transition computation, but the peak 
values were higher than either the natural or zero-length 
transition computations by more than 1597 to 20%. 
which was outside the range of experimental uncer- 
tainty. 

For r/R = 0.41 at the highest Reynolds number, 
of Re„ - 2.4.vl0 7 /m (Fig. 10), the beginning of transi- 


tion was noted downstream of the smallest penetration, 
but the boundary layer did not appear to develop into 
fully-turbulent flow. The shapes of the heating distribu- 
tions for the two larger penetrations were similar to the 
zero-length transition predictions, but the peak heating 
levels were again 1 5% to 20% higher than the computed 
values. 

Consider next the comparisons for the outer 
penetration location of r/R = 0.70 (Figs. 12-14). These 
cases are not as well suited for comparison as the r/R = 
0.41 cases, because the boundary layer did not have suf- 
ficient running length from the outer penetration for 
fully-turbulent flow to develop. Significant heating aug- 
mentation was observed in the experimental data for the 
largest hole size at all Reynolds numbers and for the two 
smaller holes at the highest Reynolds number. At Re^ = 

1.4xl0 7 /m, the heating data for the largest penetration 
fell between the natural and zero-length transition pre- 
dictions. For the two higher Reynolds numbers, the 
large penetration heating distributions appeared to 
match the shape of the zero-length transition predic- 
tions. In terms of the heating levels, the Re x = 

1.9xl0 7 /m data were approximately 10% to 15%; lower 
than the zero-length transition prediction, while the Re x 

= 2.4x1 0 7 /m data closely matched the prediction. 

The present method for predicting turbulent 
heating levels, i.e. an algebraic turbulence model with a 
specified transition location and a transition length 
bounded between zero and the running length ahead of 
the transition location, appeared to bound the experi- 
mental heating data for the range of Reynolds numbers 
and penetration sizes for the r/R = 0.70 penetration loca- 
tion. However, for the r/R = 0.41 penetration location, 
the present method under-predicted the measured peak 
heating levels by up to 20%, although the shape of the 
heating distributions was approximately reproduced. It 
should be noted that the actual heat-shield penetrations 
were not included in the computational geometry, and 
their possible effects on the flow field (i.e. circulation 
within the penetrations, shocks at the lip of the penetra- 
tions, boundary layer separation and vortex formation 
downstream of the penetrations) may be the cause for 
these differences. 

In a similar study (Ref. 5) for the Genesis mis- 
sion, the cavity was modeled in the computations. It 
was found in that study that the presence of the cavity 
had significant effects both on heating around the cavity 
itself, where a localized spike in the heating levels was 
predicted, and downstream of the cavity, where heating 
levels lower than without the cavity were predicted. 
However, it was also concluded in that study that the 
simple algebraic transition/turbulence models used (the 
same as in this study) were not sufficient to resolve the 
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effects of the cavity even with the cavity included in the 
computational geometry. 

In terms of design criteria, the worst-case for 
heating, based on the experimental data, would be a pen- 
etration at r/R = 0.41 of sufficient diameter to cause 
transition, because the boundary layer has greater run- 
ning length downstream of the penetration for turbulent 
flow to develop than it would have for a penetration at 
r/R = 0.70. This case with an r/R = 0.41 penetration 
which produced transitional/turbulent flow is repre- 
sented by the data for the two larger penetrations at Re^ 

= 2.4xl0 7 /m (Fig. 11). For this case, the turbulent 
experimental heating levels downstream of the penetra- 
tion were approximately 25% higher than at the nose. 
These experimental heating levels were approximately 
10% higher than the predicted fully-turbulent levels at 
the nose, and were approximately 20% higher than those 
predicted downstream of the penetration using the 
zero-length transition model. 


Boundary Layer Transition Correla- 
tions 


Laminar computations were performed with 
LAURA for the complete wind tunnel test range of 
angle-of- attack (a = 0-deg, 1 1-deg, 16-deg, and 20-deg) 
and Reynolds numbers (Re^ = 6.7xl() 6 /m to 2.4xl0 7 /m) 
in order to determine boundary layer edge properties 
(Re e and 8), which could be used to formulate a transi- 
tion correlation from the experimental data. Sample val- 
ues of Re g and 8 are shown in Figs. 15-17 for the Re„ = 

2.4.rl0 7 /m, a = 1 1-deg, 16-deg and 20-deg cases. 
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Figure 17: Re f) and S Distributions for Re„ = 
2.4x1 0 7 /rri, a = 20-deg 


Approximately 800 data points on the state of 
the boundary layer which covered the range of 
angles-of-attack, free stream Reynolds numbers, pene- 
tration diameters, and penetration radial locations were 
obtained during the wind tunnel test. For each data 
point, the state of the boundary layer downstream of the 
penetration was determined through visual inspection of 
the surface heating images and classified as either: lami- 
nar: localized disturbance at penetration; transition 
downstream of penetration; or fully turbulent at penetra- 
tion. Examples of experimental data which fit each of 
these classifications are shown in Fig. 18. These classi- 
fications are somewhat simplified descriptions of a com- 
plex flow field and should be used with care. For 
example, for a data point classified as “local distur- 
bance”, the heating levels measured near the penetration 
may be higher than heating levels downstream of a pen- 
etration for a data point classified as “transition down- 
stream”. Additionally, because of the three-dimensional 
nature of the flow over this geometry at angle-of-attack, 
transitional/turbulent heating levels downstream of a 
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penetration are highly dependent on the location of that 
penetration. 

In order to determine which points would be 
used to generate correlations, the data were examined to 
determine at which value of Re „ fully turbulent flow at 
the penetration was first noted for fixed values of w, r/R 
and a. To define the turbulent boundary, these points 
were fitted to an equation of the form: 

(Re v\urb = C ( >l, ' /5 )~‘ 0 ) 

where the constant was found to be C = 725. Similarly, 
the values of Re „ at which laminar flow were last noted 
(for fixed values w, r/R and a) were fitted to an equation 
of the form: 



Figure 19: Curve Fits to Wind Tunnel Data Points 
on the Boundaries of Laminar and Turbulent Flow 


= C(w/&)° ( 2 ) 

where C = 3504 and a = -2.61. The equations and the 
data to which they were fitted are shown in Fig. 19. 
Almost all of the turbulent data fell within the ±20% 
bands, but the laminar data showed more scatter. 
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Figure 18: Classification of Penetration Effect on 
Boundary Layer 



Figure 20: Entire Transition Data Set and Correla- 
tions for Laminar and Turbulent Boundaries 

Having defined correlations for the end of lam- 
inar flow and the beginning of fully turbulent flow from 
the data points on the boundaries, these correlations are 
now shown with the entire data set in Fig. 20. Although 
there is some scatter, the majority of the “laminar” and 
“fully turbulent at penetration” data points fall below or 
above the respective curves, while the “local distur- 
bance” and “transition downstream” data points fall 
between the two curves. 

In Eq. 1, the exponent on the right-hand-side 
was fixed at - 1 in order to determine a critical Reynolds 
number, as was first defined by Schiller (Ref. 22) and 
employed by many other authors (e.g. Refs. 23-26) to 
correlate transition data. Eq. I can be manipulated to 
yield: 





For wind tunnel conditions, the ratio (5/6) is 
nearly constant (~ 7.4 for this geometry). With C = 725 
and a ±20% error band, a critical Reynolds number. 
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based on boundary layer edge conditions and hole diam- 
eter, for immediate transition to fully turbulent flow at 
the penetration can then be given by: 


( Re *. t ),urb- ±20% (4) 

The boundary defined by this critical Reynolds 
number is shown for a sample case in Fig. 21. 

It should also be noted that, based on past 
experience with blunt body transition (e.g. Ref. 27), a 
conservative limit for smooth-body transition in the 
absence of penetrations has been defined as Re e > 200 
for the MSL (Ref. 4). For design purposes, it is recom- 
mended that, in addition to the smooth-body transition 
criteria, the laminar boundary correlation (Eq. 2) be 
used to determine if heating effects due to a penetration 
need to be considered for a given flight condition 
because local heating effects at the penetration may still 
be significant even if transition does not occur. 



Figure 21: Critical Reynolds Number Boundary for 
a = 16, Re„ = 9.8x1 0 6 /m, Largest Penetration Size 



x(m) 


Figure 22: Centerline Heating and /?e 0 Distributions 
for Peak Heating Point on Trajectory 
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Figure 23: Trajectory Re s and 8 Values Compared to 
Boundary-Layer Transition Data and Correlations 

Flight Peak Heating Computation 

A solution was computed for the peak heating 
point on the Mars entry trajectory (Table 2). The com- 
puted laminar heating and Re e centerline distributions 
are shown in Fig. 22. From the wind tunnel data, the 
worst-case parameters for heating were found to be the 
largest penetration diameter at the inner radial location 
(r/R = 0.4). For this case, the computed values of Re e 
and w/8 at peak heating were 263 and 5.0. For the outer 
radial location ( r/R = 0.70), the values for Re e and w/8 
were 382 and 4.3. In either case, the predicted values 
fell well above the fully-turbulent boundary from the 
correlation, and thus, turbulent heating at the penetration 
can be expected. Although not shown, the w/8 values 
for the two smaller penetrations also fell above the 
fully-turbulent boundary for both cases. Re g and w/8 
values (and heating levels) were computed at several 
additional points along the trajectory in Ref. 4. These 
values are overlaid on the transition map in Fig. 23. 
Although this correlation was generated from data 
obtained in a perfect-gas air facility and may not be 
strictly applicable to non-equilibrium, Martian atmo- 
sphere conditions, the peak heating condition was found 
to be well above the turbulent border, and thus these 
results can be taken as a conservative estimate that 
accounting must be made for penetration-induced tran- 
sition in the design of the MSL heat-shield 

Summary 

A computational and experimental study has 
been performed for the Mars Smart Lander in order to 
develop a correlation for the effects of heat shield pene- 
trations on the state of the boundary layer and to com- 
pare predicted and measured wind tunnel heating levels. 
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Heating levels also were predicted for the peak heating 
point on the Mars entry trajectory. 

Laminar Navier-Stokes predictions were found 
to agree with the laminar wind tunnel data to within the 
estimated experimental uncertainty of the data. The 
agreement in comparisons between turbulent 
Navier-Stokes predictions and experimental data was 
dependent on the location of transition in the experiment 
and the diameter of the penetration on the test model 
which produced the transition. In general, zero transi- 
tion-length and natural transition-length turbulent pre- 
dictions bounded the shape of the experimental heating 
distributions, but for the larger penetrations, the mea- 
sured heating levels exceeded the predictions by more 
than 20%. These differences were attributed to the use 
of a simple algebraic turbulence model in the computa- 
tions and the exclusion of the actual penetration (and 
hence its effects on the flow field) from the computa- 
tional geometry. 

Based on the turbulent experimental heating 
data, the worst-case heating levels were produced by the 
largest penetration size at the inner radial location: heat- 
ing levels were 25% higher than the laminar heating 
level at the nose (from experiment or computation): 
10% higher than the predicted fully-turbulent level at 
the nose: and 20% higher than the predicted zero-transi- 
tion-length heating downstream of the penetration. In 
contrast, although transition at the outer radial location 
occurred at lower free stream Reynolds numbers than 
for the inner location, there was insufficient running 
length downstream of the penetration for fully turbulent 
flow to develop and so heating levels were lower. Thus, 
for aeroheating design considerations, an outer radial 
location for a penetration would be preferable. 

Boundary layer edge quantities were extracted 
from the laminar Navier-Stokes solutions and were used 
to correlate the experimental data on the state of the 
boundary layer. The correlations were expressed in 
terms of Reg and w/5, and can be used to determine if a 
heat-shield penetration will have no effect on the bound- 
ary layer, will have a local effect, will produce transition 
downstream of the penetration, or will cause immediate 
transition to turbulent flow. Computations for the peak 
heating case on the trajectory revealed that for all pene- 
tration diameters at both the inner and outer locations on 
the leeside, fully-turbulent flow would be produced 
downstream of the penetration. 
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